function main()
close all;
h=[1,10];
k=[1,10];
epsilon=[10^-3,10^-6,10^-9];

for i=1:2
    for j=1:2
        for z=1:3
        %vettorex=GfnzNQ(h(i),k(j),epsilon(z));
        vettorex=newton(h(i),k(j),epsilon(z));
        figure,
        plot(vettorex);

        end
    end
end

end
